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We study the structural properties of self-attracting walks in d dimensions using scaling arguments 
and Monte Carlo simulations. We find evidence for a transition analogous to the transition of 
polymers. Above a critical attractive interaction u c , the walk collapses and the exponents v and k, 
characterising the scaling with time t of the mean square end-to-end distance {R 2 } ~ t 2v and the 
average number of visited sites (S) ~ t k , are universal and given by v = l/(d + l) and k = d/(d + l). 
Below u c , the walk swells and the exponents are as with no interaction, i.e. v = 1/2 for all d, k — 1/2 
for d = 1 and k — 1 for d > 2. At u c , the exponents are found to be in a different universality class. 
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In recent years different models of random walks with memory or interaction have been studied. They can be 
divided in static jl| and dynamic (^|| models, for an overview we refer to the papers of Duxbury and Queiroz Q and 
Oettinger [||. Most efforts concentrated on models with repulsive interactions, in particular self-avoiding walks (SAW), 
which have been found useful for investigating polymers in dilute solution. When an attraction term exp(— A/T), 
A < 0, is included, the SAW model reveals a swelling-collapse transition at the 'O point' T — O [0J|]. In contrast, 
the likewise challenging case of random walks with a similar attractive interaction, but without repulsion, has been 
less understood. This problem was solved only for one dimension, while in higher dimensions the results are highly 
controversial. Our numerical and analytical study of attractive random walks suggests that there exists a swelling- 
collapse transition, too, that is analogous to the transition in polymers. 

We focus on the dynamic model of self-attracting walks (SATW) 0, where a random walker jumps with probability 
p ~ exp(nu) Q to a nearest neighbor site, with n = 1 for already visited sites and n = for not visited sites. The 
interaction parameter u is equivalent to —A/T for linear polymers. For u > 0, the walk is attracted to its own 
trajectory Q|. The structural behavior of the walk can be characterized by the mean square end-to-end distance 
(R 2 (t)) and the average number of visited sites (S(t)). It is expected that these quantities scale with time t as 

(R 2 (t))~t 2 » (la) 

and 

(S(t))~t k . (lb) 



Earlier analyses for the SATW in two and three dimensions were not conclusive and the numerical data have been 
controversially interpreted p|,|s|-|lO| . While Sapozhnikov 0] considered the possibility of the existence of a critical 
attraction u c (but his numerical results were not conclusive), Lee Q and Reis Q argue strongly against the existence 
of u c , since they find v and k continuously decreasing with u. 

In this Letter we present scaling arguments and extensive numerical simulations for (R 2 ) and (S) that strongly 
suggest the existence of a critical attraction u c in d > 2, with three different universality classes for u > u c , u < u c 
and u = u c - Below u c , the SATW is in the universality class of random walks, with v = 1/2 and k = 1. Above 
u c , the SATW collapses and the exponents change to v = l/(d + 1) and k = d/(d + 1). At the critical point, the 
exponents are v c = 0.40 ± 0.01 and k c = 0.80 ± 0.01 in d = 2 and v c = 0.32 ± 0.01 and k c = 0.91 ± 0.03 in d = 3 0. 
The existence of u c is in striking similarity to the 'O point' phenomenon of linear polymers [^|,^| where three different 
universality classes for T > O, T = and T < © exist. 

We used Monte-Carlo simulations to study (R 2 (t)} and (S(t)). Figure [l] shows representative results of (R 2 (t)), for 
several values of u in d = 3. For large values of u, the curves bend down towards the slope of 2v = 0.5, while for 
small values of u, the curves bend up towards the slope of 2 v = 1. At some intermediate critical value u c = 1.9, the 
slope is approximately 1v c = 0.64. The mean number of visited sites (S(t)) shows similar behavior, with k = 1 below 
it c , k c = 0.91 at u c , and k = 0.75 above u c . Figure |^a summarizes the asymptotic exponents v and k as a function of 
u in d = 3. We obtained similar results in d — 2, the asymptotic values of v and k are presented in Fig. ||b. In Table || 
the values of the exponents are summarized and compared with the analogous known exponents for the transition 
in linear polymers. 

In the following we present analytical arguments for the exponents above criticaltity, which can explain our numerical 
findings. We assume that for sufficiently strong attraction u > u c the grown clusters are compact, so that the average 

1 /2 

number of visited sites scales with the rms displacement (R(t)) = (R 2 (t)) as 

(S{t)) ~ (R(t)) d , u » u c . (2) 

Comparing Eq. (fj) and Eq. (||) yields 

k = vd, u ^> u c . (3) 

For sufficiently strong attraction it takes a very long time for the walker to jump to an unvisited site. Before doing 
this, the walker diffuses around on the visited sites, being located with equal probability on any of the cluster sites. 
Hence the mean cluster growth rate is proportional to the ratio between the number of boundary sites and the total 
number of the cluster sites ||[l2| : 

d(S) . (R)"- 1 . t _„_ (4) 



dt 
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Thus (S) ~ t u+1 . Combining this result with Eq. ([l}d) and (||), we obtain 

1 



d+1 



(5a) 



and 



k = (5b) 



for «> u c . 

Because of universality we assume that these results, which are in agreement with the exact values v — 1/2 and 
k = 1/2 in d = 1 JTo[ ] and are supported by our extensive Monte Carlo simulations in d = 2 and d = 3, are valid for 
all u > ?i c . Indeed, Fig. || suggests that the predictions for u > u c (Eq. (||)) are approached asymptotically. We like 
to note that in d = 2 the relation k = vd also holds for u < u c , while in d = 3 the numerical results yield k < vd for 
u < u c . Since the mass of the generated clusters scales like M ~ 5 ~ R k ' v , kjv corresponds to the fractal dimension 
df of the cluster. In d = 2 the clusters are compact for all u as kjv = df = d. In d = 3 they are compact for u > u Cl 
while for u < u c , the fractal dimension of clusters generated by simple random walks df = 2 < d is obtained. At the 
criticality, we find df = 2.84 ± 0.25, but we cannot rule out the possibility that df = d. 

To understand the behavior in the critical regime we suggest the following scaling approach. Guided by Fig. |l|, we 
assume that there exists a crossover time below which the exponent v is close to v c and above which v approaches 
1/2 for u < u c and l/(d + 1) for u > u c . This suggests the following scaling relations: 

fl(t)~f*7±(t/*«) ( 6a ) 

and 

S(t) ~ t k °g ± (t/t ( ) , (6b) 

where 

t% = \u — u c \~ a . (6c) 

The plus sign refers to u > u c , the minus sign to u < u c , and the exponent a has to be determined numerically. As 
t£ is assumed to be the only relevant time scale, the scaling functions bridge the short time and the long time regime. 
To match both regimes, we require that f±(x) = const for x -C 1 (t <C t^j, and f+(x) ~ x l '^ d+1 ' Vc , f-(x) ~ x 1//2_ ' Vc 
for x ^> 1. Analogous results are expected for <7±(^), with g±(x) = const for x 1, and ,g+(a;) ~ x d ^ d+1 )~' Cc , 
9-{x) ~ x 1_fcc for x ^> 1. 

To test the scaling theory and to determine the exponent a we plotted (R 2 (t)} / and (S(t))/t^ as functions of 
i/if for several values of a in d = 2 and d = 3. We obtained the best data collapse for a — 5.0 ± 0.5 in d = 3 and 
a = 7±lind = 2, which are shown in Fig. ||a and |^b, respectively. The excellent data collapse strongly supports 
the above scaling assumptions. 
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U < U c 
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u — u c 
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SAW 

i/T = i/e 


i/T > i/e 


V 

d = 2 k 


1/2 
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0.40 ±0.01 
0.80 ±0.01 
u c = 0.88 ±0. 
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4/7 
1 

l/6»o = 0.65±0 
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1 
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0.32 ±0.01 
0.91 ±0.03 
u c = 1.92 ± 0. 
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0.59 
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1/2 
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l/6»o = 0.5 ± 0. 
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TABLE I. Comparison of the exponents v and k as well as of the estimated values for the transition points u c for random 
walks (RW) and for .SMW^ on hypercubic lattices. For values related to the O transition see M and references therein. 
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FIG. 1. The mean square end-to-end distance (R 2 (t)) versus t up to t = 10 8 timesteps averaged over 1000 configurations 
for each attraction u = 0, 1.5, 1.9, 2.25, 4 in d — 3. Note that for large values of u the curves bend down towards the slope of 
2 v = 1/2, while for small values of u the curves bend up towards the slope of 2 v — 1. 



FIG. 2. The values of the exponents k and v versus attraction u in (a) d — 3 and (b) d — 2, obtained by a least square fit 
of the slope of \n(R 2 (t)} and \n(S(t)) versus \nt for large t, respectively (see Fig. [j]). Shown are the results for t = 10 6 (A), 
t = 10 7 (v) and t = 10 s (■). Note that for u > u c and larger t the values of k and v approach the theoretical predictions of 
Eq. (^), marked as dashed lines. We estimate the value of u c to be u c = 1.92 ± 0.03 in d — 3 and u c = 0.88 ± 0.05 in d = 2 
(marked by arrows). 



FIG. 3. Scaling plots for {R 2 (t)} (O) and (S(t)) (□) for t > 1 and 20 values of < u < 3 in (a) d = 3 and (b) d = 2. For 
convenience, the data for (S(t)) have been shifted by 10 s . In d = 3 for ^ c = 0.32, k c — 0.91 and u c — 1.92 we find the best 
collapse for a = 5.0, in d — 2 for ^ c = 0.40, k c = 0.80 and u c = 0.88 we find the best collapse for a = 7. The straight lines 
represent the exponents given in Table H. 



